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ABSTRACT 

With the triaxial density profile of dark matter halos and the corresponding 
equilibrium gas distribution, we derive two-dimensional Sunyaev-Zeldovich (SZ) 
effect and X-ray surface brightness profiles for clusters of galaxies. It is found 
that the contour map of these observables can be well approximated by a series 
of concentric ellipses with scale-dependent eccentricities. The statistical distri- 
bution of their eccentricities (or equivalently axial ratios) is analyzed by taking 
into account the orientation of clusters with respect to the line of sight and the 
distribution of the axial ratios and the concentration parameters of dark mat- 
ter halos. For clusters of mass lO^^/i~^M0 at redshift z = 0, the axial ratio is 
peaked at ~ 0.9 for both SZ and X-ray profiles. For larger clusters, the devia- 
tion from circular distributions is more apparent, with rj peaked at ~ 0.85 for 
M = 10^^h~^MQ. To be more close to observations, we further study the axial- 
ratio distribution for mass-limited cluster samples with the number distribution 
of clusters at different redshifts described by a modified Press-Schechter model. 
For a mass limit of value = IO^^H'^Mq, the average axial ratio is (rj) ~ 0.84 
with a tail extended to ~ 0.6. With fast advance of high quality imaging ob- 
servations of both SZ effect and X-ray emissions, our analyses provide a useful 
way to probe cluster halo profiles and therefore to test theoretical halo-formation 
models. 

Subject headings: cosmology: theory — galaxy: cluster — large-scale structure of 
universe 



1. Introduction 



Following abundant observations of various kinds, a concordance cosmological model 
has emerged (e.g., Spergel et al. 2003). The geometry of the universe is nearly flat with 
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the total density of the universe D,tot very close to unity. Among which, about 70% is in 
the form of dark energy with negative pressure, and about 25% is composed of non-baryonic 
dark matter. The familiar baryonic matter contributes only about 5% of the total matter 
density in the universe. It is apparent that dark matter plays dominant roles in the structure 
formation. 

Besides large-scale structures, it is of great cosmological interest to understand the 
properties of individual dark matter halos. Extensive studies with high resolution simulations 
have been performed in this aspect (e.g., Fukushige & Makino 1997, 2001; Moore et al. 1998; 
Jing 2000; Jing & Suto 2000). Navarro, Frenk and White (1996, 1997) put forward a universal 
density profile (NFW profile) to describe the spherical mass distribution of dark matter halos 
of different masses. Confronted with observations, the validity of the profile has been tested 
for dwarf galaxies, galaxies, as well as clusters of galaxies. Certain degree of disagreement 
between the inferred inner profile from some observations and the NFW one raised a potential 
small-scale problem for cold dark matter models. Despite that, the NFW profile has been 
applied widely in different studies (e.g., Bartelmann et al. 2001; Li & Ostriker 2002; Wyithe 
et al. 2001). However, the sphericity assumed in the NFW profile is clearly an approximation, 
and could introduce systematic errors in theoretical predictions, and further in cosmological 
parameter extractions from observations (e.g., Birkinshaw 1999). Among others (Dubinski 
1994; Jing, et al. 1995; Thomas et al. 1998; Yoshida et al. 2000; Meneghetti et al. 2001), Jing 
and Suto (2002) did detailed analysis with numerical simulations, and proposed a NFW-hke 
triaxial density profile for dark matter halos with quantitative statistics on the distribution of 
the axial ratios and the concentration parameter. Their studies allow one to investigate the 
systematic differences caused by the non-sphericity of mass distributions. Significant effects 
have been found in the statistics of arcs, and the probability of large separation gravitational 
lenses (e.g., Oguri, Lee & Suto 2003; Meneghetti, Bartelmann & Moscardini 2003; Dalai et 
al. 2004; Oguri & Keeton 2004). 

For clusters of galaxies, the distribution of intra-cluster gas (ICM) is closely associated 
with the underlying dark matter distribution, and therefore observations, such as X-ray 
(e.g., Rosati et al. 2002) and Sunyaev-Zeldovich (SZ) effect (e.g., Carlstrom et al. 2002) 
can provide us information on the dark matter profiles (e.g., Allen et al. 2002; Grego et al. 
2004; Grego et al. 2001) . Lee & Suto (2003) studied the equihbrium gas distribution in 
triaxial dark matter halos under the approximation of small ellipticities. They (Lee & Suto 
2004) further analyzed the X-ray and SZ effect distributions associated with the triaxial mass 
distribution of clusters of galaxies. Since they aimed at finding the dark matter distribution 
for individual clusters, they parametrized the observables (X-ray and SZ) with the parameters 
being functions of those characterizing the mass profile of dark matter halos. While this type 
of reconstruction is certainly useful, individuality of clusters affects the results considerably. 
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On the other hand, for a large number of clusters, one can study collectively the profile of 
the observablcs, and further investigate the underlying dark matter distribution statistically. 
Mohr et al. (1995) analyzed 65 X-ray clusters, and presented statistical results on several 
quantities representing the morphology of clusters' X-ray emissions. With their limited 
sample, they found that the distribution of the axial ratio of X-ray emissions has a main 
peak at ?7 ~ 0.9 with a minor peak at r/ ~ 0.6. The average axial ratio is (r/) ~ 0.8. 
To investigate the three-dimensional ICM distribution, Cooray (2000) assumed a spheroidal 
ICM profile (prolate or oblate) , and obtained statistics of 3-D axial ratio given the 2-D results 
of Mohr et al. (1995). Similar studies on elliptical galaxies have been carried out for many 
years (e.g., Binney 1985; Binney & de Vaucouleurs 1981). 

The triaxial mass distribution of dark matter obtained from numerical simulations (Jing 
and Suto 2002) provides us a starting point to study the ICM distribution without artificial 
assumptions. In this paper, we analyze the corresponding surface profiles of X-ray and SZ 

effects. Instead of focusing on individual clusters, we investigate the profiles of observables 
statistically. Specifically the axial-ratio distributions of X-ray and SZ effects profiles resulting 
from the triaxial dark matter distributions are derived. Our study presents another view to 
test the theory of the formation of dark matter halos. 

The rest of the paper is organized as follows. Under the equilibrium assumption, in 
§2, wc derive the two-dimensional profiles of SZ and X-ray observables from triaxial dark 
matter distributions. In §3, we obtain the two-dimensional eccentricities of the SZ and X-ray 
profiles under the linear approximation for the deviation of sphericity. Section 4 presents 
the statistical analysis on the distribution of two-dimensional eccentricities, which can be 
compared with observational results. Summary and discussion are given in §5. Appendix A 
contains derivations of the eccentricities for the isothermal SZ effects. Results for other cases 
are also presented. Appendix B discusses possible values of the parameter C /K associated 
with X-rays and SZ effects. 

Throughout the paper, we consider the concordance cosmological model with = 0.3, 
= 0.7, Hq — 12km/ s/Mpc and (jg = 0.9, where and are the present dimensionless 
matter density and the dark energy density of the universe, respectively, Hq is the present 
Hubble constant, and cTg is rms of the extrapolated linear mass density fluctuation smoothed 
over 8h~^Mpc with h the Hubble constant in units of lOOkm/s/Mpc. 
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2. Two-dimensional SZ and X-ray profiles 

2.1. Gas density distribution 

The three-dimensional triaxial density distribution of dark matter halos is given by Jing 
& Suto (2002), which has the following form 



where R = a{x^/a'^ +1/^/6^ + z^/c^)^/^ (c < 6 < a) is the length of the major axis, Rq is 
a scale radius, 5c is a characteristic density and pcrit is the critical density of the universe. 
For the value of a, it is found that both a = 1 and a = 1.5 can fit well to the profiles of 
simulated dark matter halos. Detailed comparison showed that a = 1 is slightly better for 
cluster-scale halos while for galactic halos, a = 1.5 gives better results. In this paper, we 
focus on clusters of galaxies, and adopt a = 1.0. 

Under the equilibrium assumption, we can derive the gas distribution in a triaxial dark 
matter halo. The related equations are the Poisson equation and the equilibrium equation. 
Specifically, 

= 47rG'p , (2) 

- V = - V * , (3) 
P9 

where $ is the gravitational potential, p is the total matter density, pg is the gas density, 
and Pg is the gas pressure. 

We consider both isothermal and polytropic ICM states. For the isothermal case, we 

have 

where fc^, T^, /i, and rrip represent the Boltzmann constant, the gas temperature, the mean 
molecular weight, and the proton mass, respectively. Then the gas density is related to the 
gravitational potential through the following equation 

^ = exp[-l($-$o)], (5) 
P90 K 

where the subscript denotes the corresponding value at the central position. 
For polytropic gas, the equation of state is 

Pg = Pgo\PglpA\ (6) 
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where 7 is the polytropic index. The equihbrium gas density can be derived as 

1 7 - 1, 



Pg 



7 



(7) 



where Kq = ksTg^/ ^irip with T^o the central temperature of ICM. It is seen that the gas den- 
sity profile is closely associated with that of the gravitational potential, and the isopotential 
contours are also isodensity contours of ICM. 

As ICM contributes about 10% of total cluster mass, it is neglected in calculating the 
gravitational potential. From the triaxial density profile of dark matter halos (eq. [1]), one 
can solve the Poisson equation to obtain the gravitational potential numerically (e.g., Binney 
& Tremaine 1987). 

Under the approximation of small eccentricities for dark matter halos, Lee and Suto 
(2003) derived an analytical solution for the gravitational potential, which is given as 



$(m) ^ C[Fi{u) + F2{u) + F3(m)], 



(8) 



where u = f/Ro, f = {x,y,z) = r(sin^o cos(^o, sin ^0 sin <^o, cos^o), C — AnGScPcruRo^ , 
Cb — {I — P/a'^y/'^ and Cc — {1 — c^/a^Y^'^. Here and Cc are the two eccentricities of 
the eUipsoidal dark matter halos. The three functions Fi{u), F2{u) and ^3(14) are given in 
Lee and Suto (2003). In the following, our analytical analyses on the two-dimensional SZ 
and X-ray profiles are obtained primarily from the potential given by equation (8). On the 
other hand, we also solve for the gravitational potential numerically from Poisson equation, 
and derive the exact SZ and X-ray profiles. These exact solutions are compared with those 
analytical ones to test the validity of the approximation underlying equation (8). 

For triaxial clusters, their appearances depend on the viewing direction. With (x, y, z) 
the principal coordinate system, we introduce (a;', y', z') with z' the line of sight direction 
and x' lying in the x — y plane (e.g., Binney 1985; Oguri, Lee & Suto 2003). The two systems 
are related by 



(9) 



' x' 




' x'' 


y 


= T 


y' 


_z _ 







where 



T = 



— sin — cos cos Q cos sin ^ 
cos — sin cos Q sin sin 9 
sin Q cos Q 



(10) 
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Here 6 and are the polar coordinates of the hne of sight direction in the (x, y, z) coordinate 
system. Then we have 



C 



Fi{u') + 



-FM) 



eb'{x' cos (f) — y'sin cos 9 -\- z' sin sin + ec^(y' sin 9 -\- z' cos 



F.W) (11) 



where it' = r' / Rq. 



2.2. SZ and X-ray surface profiles 

The typical mass of a cluster of galaxies is 10^'' — 10^'^ M0, to which the gas compo- 
nent contributes about 10%. Because of the deep potential well, the typical temperature 
of ICM is a few keV . Therefore the gas is fully ionized. Hot electrons emit X-ray pri- 
marily through thermal bremsstrahlung process. On the other hand, as cosmic microwave 
background (CMB) photons pass through a cluster, they scatter with free electrons inside. 
As the result, the CMB spectrum is distorted. This is the thermal SZ effect (Sunyaev & 
Zel'dovich 1970, 1972). Observing clusters through X-ray and SZ effect has been getting 
mature, and high resolution imaging of clusters of galaxies is possible. Thus it is becom- 
ing observationally feasible to study the ICM distribution and further the underlying dark 
matter distribution in detail. 

Here we derive the respective profiles of SZ effect and X-ray emission from the gas 
distribution in a triaxial cluster. 

The CMB spectrum distortion due to the thermal SZ effect can be described by a 
frequency-dependent temperature change 5T relative to the average CMB temperature. In 
the nonrelativistic limit, the gas properties affect the SZ effect through the Compton y- 
parameter, which is proportional to the gas pressure along the line of sight. Specifically, we 
have 

5Toc JneTedl, (12) 

where rig and are the number density and temperature of electrons, respectively. The 
integration is along the hne of sight. 

For fully ionized gas primarily consisting of hydrogen and helium, Ue oc pg. Thus 



5T (X PgTedl. 



(13) 
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For isothermal ICM, the SZ effect can then be written as 



ST{x',y') oc J PgTgdl (x J — ~ ^o) d/ oc / exp 



It 

K . 



dl 



Cfe^ (x' cos (f) — y'sin (j) cos ^ + 2;' sin sin ^) ^ + (y' sin 9 -\- z' cos ^) " 



F^{u')\Yz'. (14) 
For polytropic gas, the SZ effect depends on the potential through the following equation 



(15) 



For bremsstrahlung X-ray emission, its surface brightness has a different dependence on 
the density and the temperature of electrons, which is 



(16) 



where AeH is the X-ray cooling function. For bolometric bremsstrahlung emission, we adopt 
an approximate expression AeH oc Te^^^. In the case of isothermal gas, the above equation 
is reduced to 

S,(x J nMl, (17) 

and then we get 



Sx{x',y') oc / cxp 



2C 
1< 



{x' COS (j) — y'sin (p cos 6 + z' sin (p sin 6*) ^ -|- {y' sin 6 + z' cos 9) ' 



0]}d.'. 



F3(u 

For polytropic gas, we have 

Sx{x',y')o^ \{\ 



(18) 



1 7 - 1 
7 



3+7 

[$(x',y',z')-$o]}'^'"^dy. 



(19) 



In next section, we demonstrate that both ST{x', y') and Sx{x', y') can be well described 
by a series of concentric ellipses, and we further derive an approximate analytic expression 
for their axial ratios. 
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3. Two-dimensional eccentricities of SZ effect and X-ray surface brightness 

profiles 

Lee and Suto (2003) showed that the 3-D gas distribution in a triaxial dark matter 
halo is approximately ellipsoidal with decreasing eccentricities from the center to the outer 
part of the cluster. Therefore we expect elliptical-hke profiles for both SZ effects and X-ray 
emissions. 



Let ^ be defined as: 



1 



X 



»2 



1,2 



i?o' 1 + A 1 + A" 
where the coordinates {x",y") relate to {x',y') by 



(20) 



x' 




cos ^ 


— sin \l/ 




'x"' 


.y'. 




sin ^ 


cos^ 




y. 



(21) 



The form ^ given in equation (20) represents that the length of both the major and 
minor axes is different from the circular radius f" — \J x""^ + y""^. ^Prom the view point of 
small aspherical perturbations, a circle of f" is perturbed to an ellipse labeled by ^. On the 
other hand, for an elliptical form ^ = {1/Rq)[x"'^ + y""^ /{I + A")], only one of the two axes 
is perturbed away from the circle. In other words, with same ^ and the corresponding 
original circles f" arc different. /,From Appendix A (A13), we see that ^ in equation (20) is 
in general consistent with the perturbative results. Only in the case of Cc = or 6* = 0, ^ 
gives a valid description. 

To approximate the distributions of SZ effects or X-ray emissions with a series of ellipses, 
we assume that the 2-D observables are functions of ^ only. Thus ^ represents the orienta- 
tion of an elliptical contour in {x' , y') coordinate system. With tedious but straightforward 
derivations, which are shown in Appendix A for the isothermal SZ effects, we obtain, to the 
leading order of and e^, the following general expressions for all the cases considered 



tan 2* 



Cl3 

Cn - C 



12 



(22) 



A = - 



^l+( %^^^'' )~+(Cn + C,,)\m{u'1 (23) 



-C13WI + 



A' = - 



C13WI + 



Cii — C12 



c 



13 



) +(Cii + Ci2)]m(k'), 



(24) 
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and 



1 + A 
1 + A' 



A'- 



1 + A' 




where E is the eccentricity, and 




and M is a function of u' 



+ y''^)/Ro only. The specific functional forms of M for 



different cases can be found in Appendix A. 

^Prom the above perturbative results, we see explicitly that for a relaxed cluster of 
galaxies, ^ is independent of u', i.e., the elliptical contours of an observable are in alignment 
with each other. For the eccentricity, its u' and line-of-sight {6, 0) dependences are separable. 
The -u'-dependent factor M reflects the differences between the distribution of dark matter 
and that of the ICM, and is specific for different quantities, SZ effects or X-ray emissions. 

It is clear to see that the 2-D eccentricity E is on the order of or Cc of the underlying 
dark matter halo. For a cluster, the value of E is very similar for the four considered 
cases. Figure 1 shows the eccentricity of the isothermal SZ profiles at different radius. 
The parameters of the cluster are chosen to be M = IO^^H^^Mq, 65 = 0.6, and Cc = 0.8. 
The concentration parameter is taken to be the corresponding average value (Jing & Suto 
2002). The line-of-sight direction lies in ^ = and (f) — 0. For this set of parameters, the 
perturbative result is = eb^/[—2M(u')], which is shown as the solid fine. The dotted fine 
is the exact result from direct numerical calculations without approximations. It is seen that 
the perturbative eccentricity is systematically higher than that of the numerical one, but 
not by a large fraction. The differences between the two do not exceed ~ 17% for u' up 
to 10. In terms of the axial ratio 1], the corresponding differences are less than ~ 5%. For 
smaller Cf, and Cc, the perturbative and numerical results agree better. Thus our approximate 
expressions describe the profiles of SZ effect and X-ray emission reasonably well. There is a 
mild change of E with radius. At the central region E ~ 0.5, which corresponds to the axial 
ratio r] ~ 0.87. In the region around the virial radius, E ~ 0.38 and rj ~ 0.93. 

Stark (1977) studied triaxial models of the bulge of M31. Assuming an ellipsoidal volume 
brightness distribution with constant axial ratios, x"^ + (uy)^ + {tzY = al, he derived the 
surface isophotes, and obtained a series of similar eUipses with their axial ratio depending on 
u, t and the hne-of-sight direction. Specifically for the line-of-sight direction shown in Figure 
1, the axial ratio of isophotes is (3 — 1/u, or the eccentricity a/I — (1/u)^. In the problem we 
analyze here, it is the dark matter halo that has an ellipsoidal configuration with constant 
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Cb and Cc- For the equilibrium ICM, it distributes more spherically than that of the dark 
matter with scale-dependent asphericity. Then our result E = eb^/[-2M{u')] can be clearly 
understood: Cb is the eccentricity of the dark matter halo in the line-of-sight direction, and 
the factor ^y[—2M(■u')] shows the differences between the distributions of dark matter and 
ICM. The value of E is smaller than Cb = 0.6. 



4. Statistical analysis for the eccentricity distribution 

In this section, we investigate the expected statistical distributions of E (or equivalcntly 
the axial ratio rj) of 2-D observables if dark matter halos of clusters of galaxies do have triaxial 
mass distributions as revealed by numerical simulations. Our theoretical results can be 
potentially confronted with observational results, and the agreement / disagreement between 
the two can be used as evidences supporting/challenging the cold dark matter structure 
formation scenario. 

As shown in equation (25), the eccentricity is a function of {eb,ec,Ce,0,(f)). Thus its 
statistical distribution can be calculated through the distributions of the five quantities. For 
a sample of clusters of mass M at redshift z, the differential distribution function p{E'^) has 
the form 

■2 f f f , r/dE\-^ 



dE' 



p{E^)dE^ = [- [p{c/a)d{c/a) [p{c/b \ c/a)d{c/b) I p{c,)dc, I (^) ' sinOde 
in J J J J ^ 0(p /e,ce,c/6,c/o 

(29) 

where the inner most integral reflects the projection effect, i.e., the change of the appearance 
of a fixed cluster due to different viewing. The other three integrals are for the distributions 
of the concentration parameter Ce, and the two axial ratios c/a and c/b given by Jing and 
Suto (2002). 

Before going to the full distribution of equation (29), we first study the distribution of 
E'^ due solely to the projection effect, which can be calculated through 

2 /"^^ /dE'^\ -1 

p(E^ \ce,c/b,c/a) = -/ (^7-) sinM^ (30) 

where (^i, 62) is the range of 6 in which a physically meaningful value of can be found for 
a fixed value of E'^ (Binney & de Vaucouleurs 1981). 

In Figure 2, we show p{E'^ \ Ce, c/b, c/a) of the isothermal SZ effect for a cluster of mass 
M — IO^^/i^^Mq at 2; = 0. Three sets of (e;,, Cc) are considered. In all the cases, there are two 
peaks in the distribution, which correspond to the two largest apparent axial ratios when 
viewed along the cluster's principal axes. Specifically in our definition of the coordinate 
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systems, the two line of sights are ^ = and {9 — 'K/I^cf) — 0), respectively. Our results 
are consistent with the study of Binney and dc Vaucouleurs (1981) for elliptical galaxies. 
Note that e?, and are the eccentricities of the dark matter halo and the 3-D gas profile is 
rounder than that of the dark matter distribution. 

Now let us consider the full probability function. The relevant distributions are taken 
from Jing and Suto (2002). Specifically we adopt (Oguri, Lee & Suto 2003) 

1 r |(c/a)(M„„/M.)°-°'"l')" - 0.54]^ ^ M„., ^^ "-""M" 

*'^"'=v2^xo.n3'"n amt) ■ '''' 

where is the characteristic nonlinear mass scale satisfying a'(M^) = 1.68 with cr(M^) the 
rms of linear density fluctuations smoothed over mass scale M^, and 

P(c/*|c/a) = — i— (32) 

I mm J '- ^ ^ I mm ' 

for c/fe > rmini whcrc r^nm — c/a for c/a> 0.5 and Tmin — 0.5 for c/a < 0.5. For c/h < Tmin, 
p{c/b \ c/a) — 0. For the concentration parameter Cg, we have 

/ N 1 r (Ince -Ince)^] 1 , . 

where the fltting formula of Cg in the triaxial model is given by: 



r_r 0:3 l A^irjzc) ^{z) fl + z,^V^ 

Ce - i.35exp I [(,/«)(M,,,/MO°-°^^(^)°-^J /""^V ^^^r{z) ^iz^) U + z) ^''^^ 

where Zc is the collapse redshift of the cluster of mass M^r and Amr{z) is the average density 
of a virialized halo with respect to the critical density, i.e., A^i^iz) = (SM^ir) / {4:7irl^^Pcrit) ■ 
For the numerical factor A^,, we take = 1.1. 

In Figure 3, the probability functions are shown for the four cases we discussed. To be 
more comparable with the results of Mohr et al. (1995), here and after we change p{E'^) 
to /(?]), the probabihty function of the axial ratio 77. The cluster mass is taken to be 
Myir — 10^'^h~^MQ and is placed at redshift z — 0. For all the four cases, f{r]) are very 
similar with a peak at 77 ~ 0.9 and a tail extended to 77 ~ 0.6. In the following study, we 
concentrate on f(rj) of the isothermal SZ effect. 

Figure 4 demonstrates the dependence of /(r/) on the cluster mass (left) and on the 
redshift (right). For larger clusters, the SZ profiles become more elongated, with the peak 
of f{rj) moving from 77 ~ 0.9 for M^i^ = IO^^/i'^Mq to 77 ~ 0.85 for M^ir = W^^h'^MQ. On 
the other hand, 7(77) does not show significant redshift evolution. 



- 12 - 



The above statistics focused on cluster samples with fixed masses and at fixed redshifts. 
Observationally, a cluster sample is often flux limited so that the sample contains clusters 
of various masses at different redshifts. Under the equilibrium assumption, a flux limit (SZ 
effect or X-ray) can be converted to a mass limit which is redshift dependent. Thus in 
the following we consider mass limited cluster samples. For simplicity, however, we choose 
a redshift-independent mass limit Mum for each sample. To certain extents, SZ-selected 
clusters are more or less fixed-mass-selected clusters (Holder et al. 2000). For a mass hmited 
sample, the probabihty density f{rj) can be calculated through 

. _ IdV{z)J^^^^f{v)\MAdn/dM)dM 

J{dn/dM)dMdV ' ^ ^ 

where f{i])\M.z is the probability density of t] for clusters of fixed mass M at redshift z 
analyzed above, {dn/dM)dM is the number density of clusters of mass (M, M + dM), and 
dV is the volume element. 

For the number density of clusters, we adopt the one from Jenkins et al. (2001) 

^iz,M) = 0.315^ J-I^l 
dM^ ' ' MaM dM ^ 

X exp[-|0.61 - l-a{D^aM)n (36) 

where po is the present matter density of the universe, is the linear growth factor, and ctm 
is the rms of the hnearly extrapolated-to-present matter density fiuctuation over the mass 
scale M, which is given by Viana & Liddle (1999) 

R{M) 1-7[K(M)] 



.8h-^Mpc 

where R{M) is the comoving radius corresponding to M, and 7 has the following form 



(37) 



= (0.3r + 0.2) p.92 + logio 

with r the shape parameter. 



7(i?) = (0.3r + 0.2) [2.92 + log,o (^^z^)] (38) 



In the left panel of Figure 5, the mass limited /(ry) of the isothermal SZ effect is shown for 
Miim = and = IO^^/i^^Mq, respectively. The axial ratio is calculated at the 

virial radius r^j^- It is seen that for the cluster sample with a larger Mum, the distribution 
/(?]) has a peak at a smaller 77 and extends more toward elongated configurations. The 
average axial ratio (77) f» 0.87 and 0.84 for Mum = 10^^h~^MQ and Mum = IO^^H-^Mq, 
respectively. Although only the results for the SZ effect are shown here, we expect similar 
f{r]) for X-ray emission profiles. 
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The measurements of Mohr et al.(1995) on the emission- weighted axial ratio for 57 X- 
ray clusters gave (rj) f» 0.8. The distribution has a major peak at r/ ~ 0.9 and a minor 
one at ?7 ~ 0.6. They compared the results with those of numerical simulations of different 
cosmologies. Due to the strong dynamical evolution of clusters, the profiles of X-ray emissions 
show significant non-sphericity for Q„i = 1 model with {rj) ^0.7. For low density models, the 
late time evolution of clusters is weak, and the clusters are nearly spherical with (r/) ~ 0.95 
for {Q^ = 0.2, Qa = 0) and {rj) ^ 0.91 for {Q^ = 0.2, = 0.8). The observational f{ri) 
lies in between the models of Jl^ = 1 and = 0.2. Our analysis presented here is for 
the cosmological model (Q^ = 0.3, Qa = 0.7). For M^m = lO^^h'^M^, we have (77) ^ 0.84, 
which is very close to the observational value of {rj) 0.8. On the other hand, there are 
several differences between our idealized cluster sample and the observational sample of Mohr 
et al. (1995), so that the comparison of the two is tentative. Firstly, their sample is more 
or less a fiux limited sample with some complications (Mohr et al. 1995; Edge et al. 1990), 
while our analysis is performed on a mass limited sample. Secondly, their statistics is on the 
emission- weighted axial ratio that emphasizes more on the core structure. Our results in the 
left panel of Figure 5 are for f' — r^ir. We have known that the ellipticities of the SZ effect 
and X-ray profiles are higher for the central region of a cluster than that of the outer part 
(see Figure 1). Therefore we expect smaller (r)) for the emission- weighted axial ratio in our 
model, which would be closer to the observational results. For the purpose of demonstration, 
in the right panel of Figure 5, we show the results for r' = Tmr and f' = O.Tr^r, respectively, 
for Miijn = 10^^/i~^Mq. The respective average axial ratios are (rj) 0.87 and 0.86. 

The minor peak revealed by Mohr et al. (1995) does not show up in our theoretical 
analysis. Whether this is a serious problem for the model is not clear since the number 
of observed clusters (57) is still not large enough for a solid statistical conclusion on the 
appearance of the low peak. This minor peak may just represent the tail of the distribution. 

5. Summary and Discussion 

X-ray emission and SZ effect from a cluster of galaxies are directly related to the hot 
ICM, which is gravitationally bound to the deep potential well of dark matter halos of 
clusters. Therefore, different from optically selected clusters whose identifications are affected 
considerably by the projection effects. X-ray or SZ selected clusters are relatively clean. 
Further more, the continuous distribution makes the ICM a better tracer of the overall 
underlying dark matter than galaxies. Thus X-ray and SZ observations have been used 
extensively to study the distribution of dark matter (e.g., Buote and Lewis 2004). 

Based on the triaxial dark matter halo model from numerical simulations, we derived 
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the profiles of SZ effects and X-ray emissions for clusters of galaxies. It is found that they 
can be well described by a set of concentric ellipses with the eccentricity E decreasing toward 
the outer part of a cluster. Under the approximation that the triaxiality of dark matter halos 
is weak, we obtained an analytical expression for eccentricities of these observables, which is 
a function of e?,, Cc, Cg, and 0. Here the first three quantities represent the intrinsic shape 
of dark matter halos, and the last two label the line-of-sight direction. For the four cases we 
discussed (SZ/X-ray, isothermal/polytropic), the elongation of the profiles is similar. For a 
halo of M = lO^^/i-^o and = 0.6 and = 0.8, E ~ 0.4 for ^ = and = 0. 

Our analytical results of (22) to (25) can be potentially used to constrain the dark matter 
distribution of a cluster from the corresponding observations of and E. Unfortunately, 
the line-of-sight direction (6*, 0) cannot be known beforehand, which strongly limits the 
reconstruction of dark matter distribution from X-ray and SZ observations for individual 
clusters. However, given and E, certain constraints on {9, 0) and therefore and Cc can 
be obtained. For example, as shown in Figure 1, with = 0, the observed E at the central 
region of a cluster puts a lower limit on e^. 

On the other hand, analyses on the distribution of E for a sample of clusters can 
probe the asphericity of dark matter distribution statistically. With the knowledge of the 
distributions of the five parameters, we calculate the statistics of the 2-D eccentricity or 
equivalently the axial ratio rj. We found that for clusters of mass M = IO^^H^^Mq, {rj) ~ 
0.89, and {r]) ~ 0.85 for M = lO^^h'^MQ. The distribution f{r]) depends very weakly on the 
redshift. We further studied /(?]) for mass-limited cluster samples. For M;j^ = 10^^h~^MQ, 
{t]) - 0.87. For Mum = lO^^/i^^M©, {r]) ~ 0.84. 

Keeping in mind many differences between our idealized cluster samples and the real 
observational one, we made a tentative comparison with the emission-weighted 1] distribution 
of Mohr et al. (1995) for 57 X-ray clusters. Their statistics showed that {r]) ~ 0.8. Our 
analysis for Mum = IO^^H'^Mq gave (rj) ~ 0.84 at r = Vmr- Since emission-weighted rj 
is influenced more by the central part of a cluster, a smaller (rj) is expected if a similar 
weighting is done in our theoretical analysis, which would bring a better agreement between 
the theoretical result and the observational one. 

Comparing with results of numerical simulations on different cosmologies, Mohr et al. 
(1995) concluded that low density models with — 0.2 cannot explain the observed mor- 
phology of clusters of galaxies. Although our theoretical analysis based on the concordance 
model with Jl^ = 0.3 and Qa = 0.7 shows a certain degree of consistency with that of obser- 
vations, same discrepancies as noted by Mohr et al. (1995) are apparent. The model tends 
to produce more spherical clusters than the observed ones. Gravitationally it is related to 
the weak dynamical evolution of dark matter halos for low-density models. In our study. 
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hydrostatic equilibrium is assumed for ICM, and thus the gas distribution is determined fully 
by the dark matter distribution. In reality, physical processes other than gravity can affect 
properties of ICM, and change its equilibrium profile considerably 

On the other hand, the observational sample of Mohr et al. (1995) contains 65 clusters 
(the axial ratio statistics is presented for 57 clusters), and it is still not large enough for a 
solid conclusion on the morphological statistics. Currently, a large number of high resolu- 
tion images of clusters are available. Detailed studies on their morphologies are definitely 
desirable. Because of the much better quality of images, we expect that tighter constraints 
on theoretical models of structure formation can be obtained through morphological studies 
on clusters of galaxies. 

We sincerely thank the referee for the constructive and detailed comments and sugges- 
tions. This research was supported in part by the National Science Foundation of China 
under grants 10243006 and 10373001, and by the Ministry of Science and Technology of 
China under grant TG 1999075401. 



A. Appendix 

In this appendix, we present the derivation for the 2-D eccentricity under the approxi- 
mation of small Cb and Cc- 

To the linear order of A and A', from equation (20) we have 
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where f"^ = x"'^ + y"^. In (x', y'), the above equation becomes 
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where u' = t'/Rq = y/x'"^ + y'^/-Ro, and {x',y') = (f' cos 61', f' sin 61') 



Our task is to determine the corresponding A, A' and ^ from the profiles derived from 
the triaxial dark matter distribution. For all the four cases discussed above, the derivations 
are similar. In the following, we focus on the isothermal SZ profiles. 

For 6T{x',y') = ST{^), we expand ST{^) with respect to A and A'. To the linear oder, 
we have 
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On the other hand, with f{u') = {-C / K)[Fi{u') + {l/2){eb^ + ec^)F2{u% equation (14) 
can be approximately written as: 
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The above approximation assumed that the non-circular part of ST is small due to the 
small values of and Cc- The validity of this assumption depends on the value of C/K, 
which is related to the concentration parameter of the dark halo. In Figure 6, we showed the 
comparisons between the approximate results [equation (A4) for isothermal 5T] and those 



-17- 



from the direct integrations of equations given in the previous section [equation (14) for the 
isothermal 5T] for the isothermal SZ (left) and X-rays (right) for several different values of 
C/K. It is seen that the two agree very well for reasonable C/K values, which are discussed 
in Appendix B. 

Let us define 
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Comparing with equation (A3), we obtain 
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Then to the linear order of el and e^, A, A' and ^ are determined by the following relations 
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Exphcitly, we have 
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For the other three cases, the only difference is showed up in the M factor, Mi{u') 
v! Fi2,{u') / du' Fii{u') , where i indicates different observables. 



For isothermal X-ray emission. 
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For poly tropic SZ effect. 
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For polytropic X-ray emission, 
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B. Appendix 

The parameter C/K {C / Kq for polytropic case) is related to the characteristic density 
of dark matter halos and the temperature information of ICM. Specifically, we have 

K {kBTg/jjimp) ' 

The value oi C/K depends on boundary conditions. 

Since we only consider nonspherical mass distribution up to and e^, C/K can be 
estimated under the spherical assumption. For isothermal ICM, the equilibrium equation 
tells us 

dPg _ , (r.9\ 
^ - -P'dF' ^^^^ 

and from it, we obtain 

C ^ u{d\npg)/{d\nu) 

where u = t/Rq. We adopt the boundary condition that at the virial radius, the slope of 
ICM traces that of the dark matter distribution. Then 

[1 + (2c^i,)/(l + C^ir)] 
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where C^r = Tmr/Ro. 

In the case of polytropic ICM, following Komatsu & Seljak (2001) we have 
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and 
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Fig. 1. — The eccentricity of isothermal SZ effect along the major axis. The cluster mass is 
M = lO^'^/i^^M©, and = 0.6 and = 0.8. The line-of-sight direction is ^ = and = 0. 
The solid line is the perturbative results of equation (25), and the dotted line is the results 
derived from numerical calculations without small ei, and Cc approximation. 
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Fig. 2. — The conditional probability p{E'^\ce^c/h^c/ a) for isothermal SZ effects. The mass 
of the cluster is Af = lO^'^hr^Alrij, and the redshift is 2; = 0. The eccentricity is measured at 
Tmr- Three sets of (e^, e^) are given. For each set, Cg is taken as the average value calculated 
from the fitting formula of Jing and Suto (2002). 
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Fig. 4. — The probability function for isothermal SZ effects. Left panel: The clusters 
are at z = 0. The sohd line is for M = lO^^/i^W©. The dashed line is for M = 
The dotted hue is for M = W^h'^MQ. Right panel: M = lO^'^/i-^M©. The sohd, dashed, 
and dotted lines represent results of z = 0, 0.5 and 1.0, respectively. 



-27- 




Fig. 5. — The probability function f{r}) for mass-limited cluster samples. Left panel: The 
solid hne is for Mu^ = IO^^H'^Mq, and the dotted line is for M^^ = lO^^/i-^M©. The 
axial ratio is measured at the virial radius r^r- Right panel: The solid and dotted lines 
are for the results that 77 is measured at 0.7ryir and r^r-, respectively. The cluster mass 
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Fig. 6. — The profiles of SZ effects (left panels) and X-ray emissions (right panels) along 
the major axis for isothermal ICM. The cluster parameters and the line-of- sight parameters 
are the same as those in Figure 1. The solid line in each panel represents the results from 
direct integration along the linc-of-sight [equation (14) for SZ effect]. The dotted line is 
calculated from the expansion [equation (25) for SZ effect]. The parameter C/i^ = 5, 10 and 
20, respectively. 



